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Abstract 

We present measurements of position and angular resolution of drift chambers operated 
with a Xe,C02(15%) mixture. The results are compared to Monte Carlo simulations and 
important systematic effects - in particular the dispersive nature of the absorption of transi- 
tion radiation and non-linearities - are discussed. The measurements were carried out with 
prototype drift chambers of the ALICE Transition Radiation Detector, but our findings can 
be generalized to other drift chambers with similar geometry, where the electron drift is 
perpendicular to the wire planes. 
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1 Introduction 

Around 40 years after their introduction [1], multiwire proportional chambers (MW- 
PCs) and drift chambers are widely in use in particle physics experiments and other 
fields. The main properties of these detectors, i.e. good position, timing and energy 
resolution and competitive rate capabilities at low cost, make them very attractive 
for usage in large scale high-energy physics and heavy ion experiments. 

The characteristics of these detectors have been extensively studied in the past [2]. 
However, with the stringent requirements of modern experiments and with new 
applications for proportional chambers, still a large effort is devoted to the un- 
derstanding and improvement of existing designs and to the development of new 
concepts. 

In this publication, we investigate the position reconstruction capabilities of the 
Transition Radiation Detector (TRD) [3] of the ALICE experiment. ALICE 1 is a 
dedicated heavy ion experiment to be operated at the Large Hadron Collider (LHC) 
at CERN. The ALICE TRD offers three dimensional tracking, electron/pion identi- 
fication and - combining these two capabilities - a fast trigger on high-p t electrons 
and jets. At the very high particle multiplicities anticipated in central Pb-Pb colli- 
sions (several thousand charged particles per unit of rapidity at mid-rapidity) at the 
LHC, these are very ambitious tasks. To be able to select stiff electron tracks, an 
excellent position reconstruction performance in the bending plane of the ALICE 
magnetic field is required, characterized by a position resolution below 400 fim and 
an angular resolution better than 1 °. 



2 The ALICE TRD 

In this section we describe the transition radiation detector of the ALICE experi- 
ment, in particular its position reconstruction and particle identification capabilities. 

* Corresponding author 

Email address: C . Lippmannggsi . de (C. Lippmann). 

URL: http : //www-linux . gsi . de/ "lippmann (C. Lippmann). 
1 A Large Ion Collider Experiment. 
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Fig. 1. Schematic view of a TRD chamber (not to scale). The left cross section shows a 
projection of the chamber in the x-z plane, perpendicular to the wires, the right one shows 
a projection in the x-y plane, which is the bending plane of the particles in the ALICE 
magnetic field. A particle trajectory is also sketched. The insert shows the pulse height 
versus drift time on eight cathode pads for an example event. One time bin corresponds to 
100 ns. 

2. 1 General Description and Working Principle 

The ALICE TRD consists of 540 chambers surrounding the Time Projection Cham- 
ber (TPC) in six layers at an overall length of about 7 m. The total sensitive area is 
roughly 750 m 2 ; the largest chamber is 159 cm long and 120 cm wide. Each module 
is about 13 cm thick, including radiator, electronics and cooling. The total antici- 
pated radiation thickness for six layers is about 0.15X . 

A schematic cross section of a TRD module is shown in Fig. 1. The gas volume is 
subdivided into a 3 cm drift region and a 0.7 cm amplification region, separated by 
a cathode wire grid with 0.25 cm wire pitch and 75 /im wire diameter. The anode 
wires have 0.5 cm pitch and 20 /im diameter. The drift chambers are equipped with 
cathode pads of varying sizes 2 and are read out via charge sensitive preampli- 
fiers/shapers (PASA). The whole system will consist of about 1.18 million chan- 

2 The width of the pads ranges from 0.664 to 0.8 18 cm, their length from 7.5 to 9 cm. 
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nels (readout pads). The maximum drift time is about 2 fis and the induced signal is 
sampled on all channels at 10 MHz to record the time evolution of the signal [4,5]. 
A typical signal generated by a particle track through a prototype drift chamber is 
also shown in Fig. 1. 

A 4.8 cm thick radiator is placed in front of each gas volume. This radiator is a 
sandwich of polypropylene fibers and Rohacell foam, which provides many in- 
terfaces between materials with different dielectric constants. Transition radiation 
(TR) is emitted by particles traversing the radiator with a velocity larger than a 
certain threshold [6], which for typical materials corresponds to a Lorentz factor 
of 7 1000. The produced TR photons have energies in the X-ray range (1 to 
30keV) [7] and a high-Z gas mixture (Xe, C0 2 (15 %)) is used to provide efficient 
absorption of these photons. 

2.2 Electron Identification 

The TRD will provide electron identification for momenta above 1 GeV/c [8]. To 
discriminate electrons from the large background of pions two characteristic phe- 
nomena are used: 

i) The ionization energy loss [9] at the momentum region of interest is larger for 
electrons than for pions, since here electrons are at the plateau of ionization en- 
ergy loss, while pions are minimum ionizing or on the relativistic rise. 

ii) In the momentum range considered, only electrons exceed the TR production 
threshold. 

Fig. 2 shows the mean pulse height as a function of the drift time for pions and 
electrons [4]. Here, and in the following, the time zero is arbitrarily shifted to fa- 
cilitate a simultaneous measurement of the baseline and of noise. Due to the larger 
ionization energy loss at these specific conditions (p = 2 GeV/c) the mean signal 
is about 40 % larger in the case of electrons (without radiators). With radiators the 
energy deposited by absorbed TR photons contributes considerably to the mean 
amplitude of the electrons. The characteristic signal shape for electrons with radia- 
tors is determined by the exponential probability distribution for the absorption of 
TR photons in the gas mixture. 

2.3 Tracking 

In this publication, we focus on the position reconstruction performance of the 
ALICE TRD in the bending plane of the particles in the ALICE magnetic field, 
which is parallel to the wires of the TRD and to the electric drift field. This defines 
the transverse momentum resolution of the TRD. In the third dimension, parallel 
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Fig. 2. Measured average pulse height as a function of the drift time for pions and elec- 
trons (with and without radiator). The peak at short drift times is due to the fact that elec- 
trons produced by ionization in the amplification region drift towards the anode wires from 
both sides of the wire plane, which leads to an approximate doubling of the average pulse 
height. In general, the average pulse height is larger for electrons. TR adds a significant 
energy deposit and introduces the characteristic signal shape of electrons, determined by 
the exponential absorption probability distribution of TR photons in the gas. 



to the magnetic field lines, the resolution is limited by larger pads and by the dis- 
crete wire positions. A tilted-pad design will be employed to increase the tracking 
capabilities in this direction. 

An example event in a TRD chamber is shown in Fig. 3. From the pulse height 
distribution on eight adjacent pads of 0.75 cm width, the cluster position is recon- 
structed as a function of the drift time. In this context a cluster represents electrons 
triggering avalanches in a given time bin. For a discussion of the position recon- 
struction method see section 5.3. Electrons generated close to the anode wire grid 
have a small drift time and induce signals corresponding to a small time bin num- 
ber. Electrons originating close to the drift electrode, on the other hand, have a 
larger drift time and thus correspond to larger time bin numbers. The drift time can 
be translated to a position (distance from the anode wire plane) if the drift velocity 
is known. For the ALICE TRD we aim for a drift velocity of around 1 .5 cm//xs in 
the drift region. When the clusters at all time bins are reconstructed, a straight line 
fit - as in Fig. 3 - defines the reconstructed track. The reconstructed angle (j) rec is 
obtained by: 



tan i 



aW 



V D 



(1) 
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Fig. 3. Same example event as in Fig. 1. The reconstructed clusters and a fitted track are 
overlayed. 

where a is the slope parameter of the linear fit line in pad units, W is the pad width 
in cm and is the average drift velocity of electrons in the detector in cm//is. 



3 Experimental Setup 

The measurements were carried out at momenta of 1 to 6GeV/c at the T10 se- 
condary beam line at the CERN PS. A sketch of the beam setup is shown in [10]. 
The beam was a mixture of electrons and negative pions with a momentum spread 
of about 1 %. Clean samples of each particle type were selected using coincident 
thresholds on two Cerenkov detectors and a lead-glass calorimeter. Position refe- 
rence was provided by two silicon detectors with a 50 /im strip pitch. With this 
silicon telescope the beam divergence was found to be around 0.1 (a). 

We tested four identical prototype drift chambers 3 with a construction similar to 
that for the final TRD, but with a smaller active area (25 x 32 cm 2 ). The dimensions 
of the pads were 0.75 x 8 cm 2 . We used a prototype of the PAS A with a noise on- 
detector of about 1000 electrons (r.m.s.). The FWHM of the output pulse is about 
100 ns for an input step function. The nominal gain of the PASA is 12mV/fC but 
during the present measurements we used a gain of 6mV/fC to better match to 

3 Generally, in this publication we will average over the behaviour of these four chambers, 
thus increasing the statistics of the measurements. Only where the performance of the four 
chambers is expected to be different, e.g. due to the track curvature in the magnetic field, 
we show results for a single chamber. 
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the range of the employed Flash ADC system with 0.6 V voltage swing. The high 
voltage at the anode wires was adjusted to four values corresponding to gas gains 
of 2400, 3900, 6200 and 9600. A gain value of around 4000 is anticipated as the 
nominal value for the ALICE TRD. 



4 Detector Simulations 

For simulations of the TRD performance we use AliRoot [11], the ALICE software 
package. AliRoot provides an object oriented framework for event simulations and 
reconstruction in the ALICE detector. The TRD part of AliRoot contains a full 
microscopic simulation of the detector physics processes. The interaction of the 
charged particles with the detector materials and their energy loss is simulated us- 
ing Geant 3.21 [12]. Since the production of transition radiation is not included in 
Geant 3, it was explicitly added to AliRoot. We use a momentum dependent pa- 
rameterization which applies an approximate formula for the TR yield of a regular 
stack of foils with fixed thickness, including absorption [3,13]. 

The energy transfers in primary collisions and the energy deposited by TR are con- 
verted into a number of secondary electrons and the electron collection is simu- 
lated taking into account electron drift and diffusion, amplification fluctuations, the 
distribution of the induced charge on the cathode pads (pad response), the time re- 
sponse of the detector (ion tail) as well as that of the electronics, and finally noise. 
The deposited charge is translated into raw-data-like ADC signals which then serve 
as input for track reconstruction. 

In this publication we use AliRoot to study in some detail the different contributions 
to the position reconstruction performance of the ALICE TRD. During our studies 
some changes had to be made in the AliRoot code, which will be described in the 
following. 

4. 1 Transition Radiation Absorption 

TR photons are emitted in the radiator with an angular distribution about the direc- 
tion of the emitting particle, which is sharply peaked at I/7 [6]. As a consequence, 
the TR photons cannot be separated from the incident electron track and contribute 
to the tracking information. The drift chamber detects the photoelectron ejected 
from a gas atom and the charge that is released by the secondary processes. This 
can introduce a considerable smearing of charge deposit and hence a degradation 
of the tracking performance of the detector 4 . 

4 The physical limitations imposed on the imaging quality of a xenon-filled MWPC X-ray 
imaging detector are studied in detail in [14]. 
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In the following, we only consider the xenon atoms; photon interactions with CO2 
atoms are neglected. In the absorption process, a photoelectron of energy E e = 
Ex — E s is created, where E x is the energy of the TR photon and E s is the bin- 
ding energy of the photoelectron 5 . At the relevant X-ray energies, the photoelec- 
tron is emitted preferentially in a plane perpendicular to the incoming photon track 
[15]. Even though subsequent multiple scattering and ionizing collisions with the 
gas molecules randomize the photoelectron trajectory, the charge will be deposited 
some distance away from the track. The practical range R(E e ) of this electron in a 
gas can be calculated according to [2] 



R(E e ) = AE e (l- j^TcEl) 



where A = 5.37 • 1(T 4 gcm- 2 keV-\ B = 0.9815 and C = 3.123 • lO^keV -1 . 
As an example, for our gas mixture the range of a lOkeV electron is about 500 fim. 

The emission of the photoelectron leaves a hole in the shell which will be filled 
with electrons from higher shells. This de-excitation occurs by emission of either 
an Auger electron or a fluorescence photon. The probability for de-excitation by 
photon emission is determined by the fluorescence yield, which is 0.87 for the 
xenon K-shell [16]. The fluorescence photon energy is Es — E T , where E T is the 
binding energy of the second shell involved in the transition. The emission of fluo- 
rescence photons is isotropic and their absorption length in the gas is exponentially 
distributed with a mean that is given by the attenuation coefficient. As an example, 
a K-shell fluorescence photon will have an energy of E K 35keV; in pure xenon 
it will have an absorption length of 24 cm and can generate a background hit at a 
distant position. However, since the energy of the largest part of the TR photons is 
below E K , these are rare events (i.e. 1.4% at p = 2GeV/c). L-shell fluorescence 
X-rays, on the other hand, are very common (due to the mean TR photon energies 
around lOkeV and the large fluorescence yield); they carry an energy of around 
5 keV and have an absorption length around 0.4 cm in Xenon. Auger electron emis- 
sion is also isotropic. The range of the Auger electrons is calculated by Eq. 2. As an 
example, an Auger electron emerging from the L-shell (E K — 2E L 25 keV) has a 
range around 0.2 cm in Xenon. A simplified picture of the just described secondary 
processes has been added to the AliRoot code to allow studies of their influence on 
the position reconstruction performance of the TRD. 



5 The xenon K-shell binding energy is about 35 keV, the average L- and M-shell binding 
energies are about 5.1 and 0.9 keV, respectively. 
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Fig. 4. (a) Average drift time tjj for electrons at different positions in the drift volume. The 
anode wire is situated at (x = cm, z' = cm), z' is the distance to the nearest anode wire 
in z (Figs. 1,5). The calculation was done with GARFIELD for the Xe,C0 2 (15 %) gas 
mixture, an anode wire voltage of 1550 V and a drift voltage of -1950 V. (b) Corresponding 
distribution of the average drift velocity. Note that the coordinate system has been rotated 
as compared to (a) for better visibility. 

4. 2 Electron Drift Path 

In drift chambers one takes advantage of a unique relation between the position of 
primary ionization electrons (x) and the drift time (to) to the nearest anode wire, 
where the electrons generate avalanches. To precisely reconstruct the position of 
the passage of the particle through the detector one generally wants to know this 
space-time relation, which may not be linear: 




(3) 



o 



Here vo(t) is the local drift velocity at time t. For constant drift velocity the space- 
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Fig. 5. Ideal drift lines for electrons coming from the drift region (from the top) at different 
z positions. The anode wire plane is at x = cm; the cathode wire plane is at x = 0.35 cm. 

time relation Eq. 3 becomes linear: 



x = v D t D . (4) 

In the TRD v d is constant in a large fraction of the detector (vd — vq = 1.5 cm//is 
in the drift region), but (in general) higher in the amplification region. However, 
one can approximate: 



x « v$t D , (5) 

where vf? is an average drift velocity. We used GARFIELD [17] to calculate drift 
times to (Fig. 4a) and average drift velocities (Fig. 4b) for electrons generated 
at a given (x, z') position. Here z' denotes the lateral distance of the position of the 
drifting electrons to the closest anode wire (0 < z' < 0.25 cm). A large number of 
electrons were drifted from each position (including diffusion) and the most proba- 
ble drift time was taken. We find that vf? is approximately equal to v only close to 
the drift cathode (x ~ 3.35 cm), but increases as x decreases. In the amplification 
region it is in general more than twice as large (3.5 cm//is). 

However, the drift time to depends also on z' . Fig. 5 shows the different drift paths 
for electrons coming from the drift region at different ^-positions. At a given dis- 
tance from the anode wire plane (which is situated at x = cm), the shortest drift 
time is given by electrons drifting at z' = cm. Electrons drifting at z' = 0.25 cm 
have a longer drift path and - on top of that - cross the low field region between 
two anode wires. Consequently, we observe in Fig. 4a a drift time offset depending 
on z' . For electrons coming from the drift cathode (x ~ 3.35 cm) it is around 120 ns 
as compared to the value at (x = cm, z' ~ cm). For electrons originating in the 
low field region at (x = cm, z' = 0.25 cm) we find v™ < 1 cm//is, resulting in a 
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drift time offset 6 up to 430 ns! 



The non-linearity of the space-time relationship as illustrated in Fig. 4 has been 
added to the AliRoot code. 



4.3 Pad Response Function 

Proportional chambers often feature a cathode plane subdivided into separate strips 
or - like the ALICE TRD - pads with independent charge sensitive readout for the 
purpose of localizing the avalanche with a precision that is a fraction of the strip or 
pad width W. A parameter that strongly influences the distribution of the induced 
charge on the cathode plane is the angular position of the avalanche at the anode 
wire [18]. However, in most practical situations it is not possible or desirable to 
restrict or control this quantity so that one generally observes a resultant effect due 
to all avalanche angles. There exists an empirical formula for the induced charge 
distribution p(y) by Mathieson [19] that describes well such average behaviour in 
symmetric MWPCs along the anode wires 7 . The coordinate y is given by the wire 
direction (Fig. 1). 

The pad response function (PRF) P(y) is obtained by integration of p(y) over the 
width of the strip or pad 



The ALICE TRD however is not exactly a MWPC, but is extended by a drift volume 
that is separated from the amplification volume by a cathode wire grid. As we shall 
see in this section, the Mathieson formula can nevertheless be used to calculate 
P(y) for this design to a rather good accuracy. 

The exact PRF P(y) can be obtained by employing the weighting field formalism. 

The weighting field Ew(y) is the (imagined) electric field in the detector when 
the readout electrode is set to 1 V while all other electrodes are grounded. The 
field E w (y) is generally used to calculate induced currents in arbitrary electrode 



Diffusion considerably broadens the drift time distributions, especially in this region. 
Fig. 4 only shows the most probable drift times. 

7 A symmetric MWPC consists of a plane of anode wires centered between two planar 
cathodes. 



y+W/2 




(6) 



y-W/2 
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Fig. 6. Pad response functions for W = 0.75 cm. The circles show the measured PRF, the 
dashed line shows the results of a calculation using the Mathieson formula and the solid 
line shows the exact calculation for the ALICE TRD geometry. 

geometries, using the Ramo theorem 8 [20] via 

i{t) = -qE w (-r>(t))-l?(t) . (7) 

The current i that is induced on a readout electrode at time t by a charge q moving 
with velocity ~v is proportional to the weighting field E w at the position of the 
charge. Calculating E w as a function of y on the cathode plane yields the cathode 
charge distribution p(y) for a given geometry. From this we can then calculate P(y) 
following Eq. 6. 

We used GARFIELD to calculate p(y) for the ALICE TRD geometry (Fig. 6). 
Since we are - as already mentioned - not interested in the influence of the angular 
position of the avalanche around the anode wire, we average over many angles. 

The exact method confirms that the Mathieson formula is a good approximation. 
In general, the PRF can be well approximated by a Gaussian curve [2]. The stan- 
dard deviations op of Gaussian fits to the P(y) (in pad units) are 0.495 and 0.496 
for the measured PRF and the exact calculation for the ALICE TRD geometry, re- 
spectively. For the Mathieson formula and for the exact calculation for a MWPC 
geometry with similar parameters 9 we find 0.482 and 0.485, respectively. The PRF 
for a symmetric MWPC is narrower by about 3 % (in a P ) as compared to the TRD 
geometry, for our specific wire diameters, wire pitches and anode-cathode separa- 
tion. The PRFs calculated with GARFIELD, as described in this section, are used 
for the simulation of the pad response in AliRoot. 

8 Also known as the reciprocity theorem. 

9 No drift region. 
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5 Position Reconstruction and Systematic Effects 

In this section we present some first results on the position reconstruction perfor- 
mance of the ALICE TRD and describe the different systematic effects. 

5. 1 Definitions of Resolutions 

The residuals for a given track are defined as the distance between the position of 
the reconstructed cluster (y t )d and the position of the reconstructed track (y t ) f it for 
each time bin t: 

A y = (y t )d ~ {Vtjfit ■ (8) 

As the position resolution a y of the detector we define the sigma of a Gaussian fit 
(within 3a) to the distribution of residuals A y for a large number of tracks. This 
resolution does not depend on 'external effects' like multiple scattering in front of 
the gas volume of the drift chamber and/or beam divergence. It thus represents the 
detector-intrinsic position resolution. As the angular resolution we define the width 
00 of a Gaussian fit (within 3a) to the distributions of the reconstructed angles. 
This resolution includes the mentioned external effects. 

5.2 Tail Cancellation 

The signals that are read out from the cathode pads are induced by the positive 
ions generated in the electron avalanches near the anode wires. Since the massive 
ions move slowly compared to the electrons, the signals exhibit long tails. Con- 
volution with the response of the PASA yields the time response function (TRF), 
which is asymmetric. For our specific chamber geometry and electronics, the tail 
of the TRF can be well approximated by the sum of two exponential functions with 
characteristic decay times 

T short « 0.10 /is and T long « 0.93 /is . (9) 

The TRF gives rise to a strong correlation between the signal amplitude in sub- 
sequent time bins. This is in general a problem also in other related detectors, in 
particular in TPCs, since it biases the position measurement results as a function of 
time. In the case of the TRD the correlations affect especially the angle measure- 
ment (see Fig. 3). A way to minimize the effect is to remove the tails from the data 
by deconvolution (tail cancellation). Three different methods are studied here: 
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Fig. 7. Time dependence of the average PASA pulse height for 55 Fe X-rays. Due to the 
pointlike charge deposit by X-rays this signal is almost identical to the time response func- 
tion (TRF). In the shown signal diffusion is included, in the TRF not. The upper panel 
shows the original signal and the signal after applying a tail cancellation with one and two 
exponential functions (ExplTC, Exp2TC) . The lower panel shows the effect of adding the 
short tail component to the left (tail making, TM), and of subsequent tail cancellation (long 
component). 

• The one exponential tail cancellation (ExplTC) subtracts the tail (for each time 
bin) as a function of time. Here the tail is assumed to be a one exponential func- 
tion with decay time Ti ong . 

• The two exponential tail cancellation (Exp2TC) subtracts accordingly a tail that 
is assumed to consist of a superposition two exponential functions with decay 
times Ti ong and T short . 

• Finally we also apply a signal symmetrization (TM+TC) that first replicates the 
tail with T s h ort at the times preceding the maximum (tail maker, TM) and then 
subtracts only the long component similar to the first mentioned method (TC). 

The effects of the three different methods on the TRF are shown in Fig. 7. The 
tail of the original TRF is largely reduced by the ExplTC method and the maxi- 
mum signal amplitude is lowered by around 10%. However, the TRF is not fully 
symmetrized, so we expect some correlation to remain, if this method is used. The 
Exp2TC method symmetrizes the TRF but the effective signal amplitude is reduced 
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Fig. 8. Measured histograms of residuals A y (left panels) and of the reconstructed an- 
gles (right panels) for 6 GeV/c pions. We show data without tail cancellation and for the 
ExplTC, Exp2TC and TM+TC methods (from top to bottom). Where applicable, Gaussian 
fits are also shown. 

by around 30 %, which introduces a considerable degradation in the signal-to-noise 
ratio. The TM+TC method also symmetrizes the TRF but without the drawback of 
a reduction in signal amplitude. 

Fig. 8 shows example histograms of the residuals A y of fitted tracks and of the 
reconstructed angles for the different tail cancellation methods described in sec- 
tion 5.2. The incident angle of the beam was ~ 15°. If no tail cancellation is 
applied, the distributions are very broad, with pronounced tails as a result of the 
mentioned correlations. Since an entry at a given time bin increases the amplitude 
at later time bins on the same pad through the TRF, the reconstructed angles are 
generally shifted towards smaller values by this effect. Tail cancellation improves 
the situation but in the case of the ExplTC method the correlation is not fully re- 
moved. However, the correlations are nicely removed by the Exp2TC and TM+TC 
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methods. 



The best results at the described conditions are obtained using the TM+TC method 
(o-y = 232 /im and = 0.687 °). With the ExplTC method we find a y = 337 fj,m 
and = 0.686°. 

For Nfu independent fitpoints the relation between the accuracy of the measure- 
ment in y of the single points (here represented by a y ) and of the angular resolution 
&<t> is given by [2] 



/ 12 (Xy 

a * *> iw* d rad ■ (10) 

Following Eq. 10, we expect for about Nf it = 20 fitpoints (see Fig. 3), for a detector 
thickness of 3.7 cm and for a position resolution of a y = 337 yum - as measured 
with the ExplTC method - an angular resolution of o$ = 0.4°. For the TM+TC 
method we expect an even better angular resolution: = 0.23°. The measured 
angular resolution does not reach these expected numbers. We conclude that the 
fitpoints are not independent as presumed by Eq. 10, since the different correction 
methods supposedly do not remove fully the correlations between the signals in 
subsequent time bins (between the fitpoints). 



5.3 Cluster Reconstruction 



For each time bin charge sharing between adjacent pads allows to reconstruct the 
position of the clusters along a pad row (in the wire direction). To calculate the y- 
position of that cluster (Fig. 1) we assume a Gaussian PRF 10 . The amplitudes in at 
least two neighbouring pads are required to be above threshold, which is determined 
by the value of the noise N. N was extracted from the baseline in the presamples 
of the drift chamber signals. Gaussian fits to the noise distributions yield values 
of around 1.7 ADC channels. The displacement y^ is of the cluster from pad % is 
calculated using a weighted mean of two measurements [2]: 



1 

Vdis 



Wi + U> 2 



<T I , Ai 



wi I ^ In 



W M-1 




. (11) 



Here ap is the Gaussian width of the PRF, W is the pad width and w\, W2 are 
weights: wi = (Aj_i) 2 , w 2 = (A i+ i) 2 , with Ai being the amplitude on pad i. The 



10 A simpler method is to calculate the center of gravity on three pads, but the PRF method 
yields results that are more accurate by about 10 %. 
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error of the cluster position is given by 



a t = 50 ns in time direction and (12a) 



v y = \/Vo) 2 + -t^ cm along?/. (12b) 



Here A = + Ai + A i+ \ is the sum of the amplitudes on the three pads, with 
Ai > Ai-\ an d Ai > A i+1 . The parameter a rs 0.03 cm is a specific resolution 
that is optimized for the best detector performance. It is of the order of the residuals 
A y . We also apply a center of gravity correction to the time coordinate (x in Fig. 
1), by moving the reconstructed clusters in that coordinate according to the values 
of the amplitudes in the neighbouring time bins. Assuming on a given pad the three 
amplitudes A t -i, A t and A t+ ± at three subsequent time bins, the cluster at time bin 
number t is shifted to 



t + St = t + 5 ~ At -\ +At ": 1 • (13) 
At-i + A t + A t+1 

5o ~ 1.5 is a factor that is optimized for best detector resolution 11 . This proce- 
dure corrects for the ambiguity of the position of the signal within the 100 ns time 
bins and significantly improves the position resolution: From Fig. 9 we find an im- 
provement of 35 %. Despite the better position resolution, the angular resolution 
is improved by only 1 %. Following Eq. 10 this would again imply that to some 
extend the correlations between the signals in subsequent time bins cannot be re- 
moved. The obtained angular resolutions seem to be the lowest limit. 



5.4 Non-Linearities 



Fig. 10 shows the systematic variation of the reconstructed angle rec with the z- 
coordinate across the wires, extracted from the silicon strip detectors. Clearly visi- 
ble is the influence of the anode wire grid with 0.5 cm periodicity. This systematic 
effect can be approximately reproduced by the simulation and is explained by the 
non-linearities in the time-space relationship as discussed in section 4.2. A small 
variation of the angle 9 from zero has to be assumed (here 9 = 1.5 °), indicating a 
slight misalignment of the chamber with respect to the beam. In fact, the precision 

11 Since now the distance in time direction between reconstructed clusters is not constant 
anymore, it seems necessary to change the value of the error at given by Eq. 12a. However, 
the effect on the resolution obtained is small and we keep at constant, as in Eq. 12a. The 
result of the time center of gravity correction can be seen in Fig. 3. 
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Fig. 9. Measured histograms of residuals A y (top panel) and of reconstructed angles (lower 
panel). We compare data with (solid lines) and without (broken lines) time center of gravity 
correction (COG corr.). 



of the alignment in this direction was of this order. If this is the case, the lateral 
distance of electrons deposited along the tracks from the nearest anode wire in the 
drift cell (z' in Fig. 4) is varying with the distances from the anode wire plane x. 
Thus an offset that depends on x is added to the drift time of the electrons, introduc- 
ing the observed systematic effect. The resolution deterioration due to this effect is 
about 0.36 °, at these specific conditions (pions, 3 GeV/c, w 15 °). For 9 = 0° the 
systematic effect disappears in the simulated data. 



Fig. 1 1 shows the residuals A y as a function of the position of the avalanche with re- 
spect to the pad. This value is a measure of the error of the coordinate measurement 
using Eq. 1 1 . For avalanches in the center of a pad the error in the measurement is 
about 20 % larger. Since the same reconstruction method is used in the simulations, 
this effect is also well reproduced. 
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Fig. 10. Systematic variation of the reconstructed angle <p rec with the z-coordinate (across 
the wires). We show measured (crosses) and simulated results (solid line). The insert 
shows a projection of the measured data on the ordinate, giving the overall effect due to 
non-linearity. The TM+TC method was used. 
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Fig. 11. Residuals A y of the fit as a function of the position of the avalanche with respect 
to the pad. We show measured (crosses) and simulated results (solid line). The TM+TC 
method was used. W is the pad width. 



6 Position and Angular Resolution 



Here we present the measured detector performance and compare it to AliRoot 
simulations. 



19 



0.07^ 
0.06r 
'E 0.05- 
0.04r 
0.03r 
0.02r 




H — | — I — I — I — | — I — H 



1; 

„ 0.9: 

CD 

-g 0.8; 







0.7; 
0.6 : 
0.5: 



t 1 1 1 1 1 1 1 1 1 1 1 1 r- 

■ e, Data 

• 71, Data 

□ e, Data, no Radiators 

O jc, Data, no Radiators 

e, Simulation 

it, Simulation 

e, Simulation, no Rad.- 



H — I — I — I — I — I — I — I — I — I — H 

p = 3 GeV/c 
<t> = 1 5.5 deg. 
6 sim =1.0deg. " 




20 40 



60 80 
S/N 



100 120 140 



Fig. 12. Position resolution a y and angular resolution as a function of the signal-to-noise 
ratio (S/N). The filled squares (circles) show the measured data for four different gain 
values for electrons (pions) with radiators in front of the drift chambers. The large open 
squares (circles) show the measured data for electrons (pions) with no radiators and for one 
given gain. The lines show simulation results for different gains for pions and electrons 
with and without radiators. 



6.1 Dependence on S / N 



In this section we study the performance of the detector as a function of the signal- 
to-noise ratio (S/N). The signal height S was extracted from pulse height spectra 
(mean value) at a given time bin corresponding to the center of the drift region 
(1.6 /is drift time at nominal conditions as in Fig. 2). 

The dependence of a y and 0$ on S/N is shown in Fig. 12. Again the incident angle 
of the beam was about <p = 15 °. The measured data are nicely reproduced by the 
AliRoot simulations for 9 = 1 (see Fig. 1). The measured data points for pions 
and electrons lie on two separate curves, roughly of 1/ y 1 S/N form. At a given S/N 
value, the resolution is worse for electrons as compared to pions. Since the S/N 
value at a given gas gain is about 60% larger for electrons, at normal operation con- 
ditions the resolution is very similar for both particle types. The data points without 
radiators shows better resolutions for electrons and lies on the same curve as the 
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pion data , while for pions the performance is similar with and without radiators. 
This indicates that the deterioration of the resolution in the case of electrons with 
radiators is connected with one of the following two processes: 

i) Bremsstrahlung created in the radiator, 

ii) Transition Radiation from the radiator. 

The simulations reproduce the observed behaviour well, implying that the brems- 
strahlung contribution is very small. The processes secondary to the TR absorp- 
tion on the other hand turn out to cause a significant deterioration of the detector 
resolution for the electrons. As described in section 4.1, L- shell fluorescence X- 
rays are very common. They carry an energy of about 5 keV and their mean free 
path is about 0.4 cm in Xenon. Their influence on the resolution is dominant; the 
influence of the range of photoelectrons and Auger electrons on the other hand is 
small 13 , as well as the influence of K- and M-shell fluorescence X-rays. The former 
are high-energetic and generally escape from the region where the TR absorption 
takes place. The latter are low-energetic and their absorption length is too small 
100 /im) to effectively influence the resolution. 



6.2 Dependence on Incident Angle 

Fig. 13 shows the position resolution a y and the angular resolution as a function 
of the reconstructed angle (p rec for two tail cancellation methods. The measured re- 
sults are quantitatively reproduced by the AliRoot simulations. Only at small angles 
is the simulated resolution generally better than the measurements. This systematic 
deviation can be explained by a space charge effect [21]. At small values of the 
incident angle (perpendicular tracks) all electrons created along a track drift to the 
same anode wire spot 14 , leading to a buildup of positive ions around this wire spot. 
This reduces the effective gain and thus the S/N value for these conditions. As a 
consequence, we observe a deterioration of the resolution for small angles, which 
is not reproduced by the simulations, since no space charge effect is included. 

For the ExplTC method the measured position resolution is around 350 /im at 
4> rec ~ 15 ° and improves for smaller values of the incident angle. For the TM+TC 
method the position resolution is below 250 fim at all investigated angles. The angu- 
lar resolution is around 0.7 ° at <p rec ps 15 ° and below that value for smaller incident 

12 The S/N value for electrons decreases without radiators due to the absence of energy 
deposit by TR. 

13 This agrees well with our observation that a magnetic field of 0. 14 to 0.56 T has no 
influence on the resolution. The tracks of photo and Auger electrons would be curled up in 
magnetic fields of that strength, which would lead to an improved resolution. 

14 If a magnetic field is applied, the angle <f> where the space charge effect is largest is 
modified by the Lorentz angle 4>L (see Eq. 14). 
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Fig. 13. Position resolution a y (upper panels) and angular resolution o$ (lower panels) as a 
function of the reconstructed angle for the ExplTC (left panels) and TM+TC (right panels) 
methods. The radiators were mounted in front of the drift chambers and the simulated data 
correspondingly contains TR in the case of electrons. 



angles and both tail cancellation methods. While the TM+TC method does not im- 
prove the measured angular resolution - as was already mentioned - we observe a 
slight improvement in the simulated angular resolution. The observed performance 
is well within the requirements for the ALICE TRD that were listed in section 1. 



6.3 Dependence on Drift Velocity 



The nominal drift velocity for the ALICE TRD of around vq = 1.5cm///s in the 
drift region was determined together with the sampling rate of 10 MHz to lead to a 
sufficiently large number of fit points for the track reconstruction (Nf it ~ 20, see 
Fig. 3). However, to provide more general results we also varied v D and investigated 
its influence on the resolution of the detector. Fig. 14 shows the dependence of the 
position resolution a y and of the angular resolution on v^,. In the measurements 
the average drift velocity ~ vrj can be extracted from the data using Eq. 1, if the 
incident angle of the beam is known. Values for vd were varied by setting the drift 
voltage to 2.1 (nominal value), 2.4 and 2.7 kV, while keeping the anode wire voltage 
constant at 1 .55 kV, corresponding to a gain of 3900. The two main contributions to 
the resulting resolution values are the lever arm of the fit (number of fit points) and 
the average amplitude per time bin. The number of fit points depends inversely on 
the drift velocity. The average amplitude per time bin is increased by a larger drift 
velocity, since more electrons reach the anode wires per time unit, leading to an 
increase of the S/N value. As a consequence we find, in Fig. 14, that a large drift 
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Fig. 14. Measured and simulated dependence of the detector performance on the drift ve- 
locity for electrons and pions. From top to bottom the position resolution a y , angular reso- 
lution cr^ and signal-to-noise ratio S/N are shown. The drift velocity scan was carried out 
with radiators. 



velocity leads to a deterioration in the resolution due to the decrease in the number 
of fit points. Accordingly, the resolution also deteriorates for small drift velocities, 
due to the decrease of S/N. Our nominal conditions in the beam test - namely a 
drift voltage of 2.1 kV, corresponding to a drift velocity of around v D = 1.67cm//is 
- turn out to be a good choice for our specific chamber dimensions and readout rate 
(10MHz). 



6.4 Dependence on Momentum 



In Fig. 15 we show the dependence of the resolutions o y and on the beam mo- 
mentum p. Generally, o y and for pions improve for larger momenta, which is 
explained by the increased S/N value for larger momenta. The ratio S/N as a func- 
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Fig. 15. Measured and simulated dependence of the detector performance on the particle 
momentum for electrons and pions. The position resolution a y , angular resolution and 
signal-to-noise ratio S/N are shown (from top to bottom). The experimental momentum 
scan was earned out without radiators. However, here we added also simulated data with 
TR to demonstrate its momentum-dependent effect on the resolution of the electrons. 



tion of p is shown in the lower panel of Fig. 15. At p w 0.56 GeV/c, corresponding 
to j3j = p/(m n c) m 4 (where m n ~ 140 MeV is the mass of the pion), pions 
are minimum ionizing; at larger momenta - especially in the momentum range of 
interest between 1 and 6 GeV/c - the ionization energy loss and thus the measured 
S/N value continuously increase [9]. Electrons already have a factor /3j m 2000 
atp = 1 GeV/c momentum. Thus they are at the plateau of ionization energy loss. 
In the momentum range of interest, the S/N value is constant (with no radiators), 
leading to a constant resolution as a function of the momentum. The situation is 
changed with radiators, since at an electron momentum of p = 1 GeV/c the TR 
production sets in and at higher momenta the energy deposit due to TR and thus the 
S/N value increase considerably. However, due to the effects described in section 
6.2 the performance is not improved by the larger energy deposit associated with 
TR, but deteriorated. 
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6.5 Performance in Magnetic Field 



The ALICE TRD will be situated inside the large L3 barrel magnet with a mag- 
netic field of 0.4 T. Since the electrons drift perpendicularly to the field, they will 
experience the Lorentz force that leads to a displacement of the clusters along the 
pad rows as a function of the drift time. For an electron produced at position (x , 
Do, z ), where the y and z directions are perpendicular to the drift direction along x 
and the y direction is parallel to the wires (see Fig. 1), the new y position is given 
by 



y = y + lut(x-x ), lut = tan0 L , (14) 

where 4> L is the Lorentz angle. It is visible as an apparent inclination of the recon- 
structed track. If the track passing through the detector has an angle 4> , then the 
reconstructed angle <p rec is given by 



4>rec = 4>0 ~ 4>L + $4> , (15) 

where 8(f) is the error of the measurement. The Lorentz angle depends on the mag- 
netic field strength and the drift velocity of the electrons. This dependence needs to 
be known to be able to reconstruct the original particle track and extract its inclina- 
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Fig. 17. Position resolution a y (top panel) and angular resolution (lower panel) as a 
function of the reconstructed angle with and without magnetic field for pions. The symbols 
show the measurements, the lines show simulated data. In the data with magnetic field the 
actual incident beam angle was around 4>q = 2.2° . The magnetic field dependent Lorentz 
angle is added to <fi following Eq. 15. We also show measured data taken with a 0.8 cm 
thick aluminum plate in front of the detectors (stars and crosses). 



tion 0o- 

We measured the Lorentz angle as a function of the magnetic field and of the elec- 
tron drift velocity and compare the results to MAGBOLTZ [22] calculations (Fig. 
16). For a magnetic field of 0.4 T and a drift field of 750 V/cm we find a Lorentz 
angle close to 8 °. 

In Fig. 17 we show the dependence of the resolution on the Lorentz angle <pL. It 
is very similar to the dependence on <p , which implies that the resolution depends 
only on the value of reconstructed angle rec , and that additional effects due to 
the presence of the magnetic field are negligible. In Fig. 17 we also show data 
taken with a 0.8 cm thick aluminum plate in front of the detectors. This plate has 
a radiation length of about 0.9X , similar to that of four TRD layers. As expected, 
the resolution in a y is undisturbed, since it does not depend on external effects like 
multiple scattering. On the other hand, we find clear effects on evidencing a 
momentum dependent multiple scattering. 
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6. 6 Performance compared to External Track Reference 
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Fig. 18. Position resolution of the detectors with respect to a silicon strip detector for the 
ExplTC (upper panel) and TM+TC (lower panel) methods for pions at 4GeV/c. The ex- 
tracted position resolution of 376 [im for both methods is for the reconstructed tracks at the 
center of the detector (time bin 15, 1.5 ^s drift time, see Fig. 3). It includes the resolution 
of the silicon strip telescope (« 50 ^m) as well as the beam divergence (« 0.1 °). 



Finally we show in Fig. 18 the position resolution of the detector with respect to 
a silicon strip detector, which is located a few centimeters in front of the investi- 
gated chamber. Let y^c be the position of the center of the reconstructed track (at 
time bin 15, see Fig. 3) and y si the position given by the silicon telescope. Then 
Vsi — Vdc is a measure of the total position resolution of the drift chamber. This 
does however include the resolution of the silicon strip detector (?» 50 yum) as well 
as the beam divergence (~ 0.1 °) and external scattering effects. Fig. 18 shows a 
histogram of y si — y DC for pions at p = 4 GeV/c, at O « 15 °, and for two tail 
cancellation methods. A Gaussian fit yields a position resolutions of 376 /im for 
both methods, which is similar to the resolution extracted from the residuals for the 
same conditions. 
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7 Summary and Conclusions 



We have measured the performance of drift chamber prototypes for the ALICE 
Transition Radiation Detector (TRD) with respect to position and angular reso- 
lution. The detectors are drift chambers with cathode pad readout filled with the 
Xe,C02(15%) mixture. For incident particle angles from to 15 ° with respect to 
the wire normal we find a position resolution better than 300 /im (a) and an angular 
resolution below 0.8° (a). A systematic effect of about 0.36° at <p = 15° is in- 
troduced by non-linearities: The discrete configuration of the wire grids in connec- 
tion with the generally higher drift velocity in the amplification region introduces 
a modulation in the electron drift times, leading to a distortion of the space-time 
relation (non-linearity). 

If a radiator is added to the drift chambers, transition radiation contributes to the 
energy deposit in the gas. Then the signal-to-noise ratio is increased for electrons, 
but, nevertheless, the electron resolution is by about 7 % worse in that case. L- 
shell fluorescence photons, which are produced in secondary processes after the 
absorption of the transition radiation photons, have an absorption length of about 
0.4 cm in the xenon gas mixture. This smearing of the charge deposit around the 
actual track of the incident electrons introduces a considerable degradation of the 
position reconstruction performance for electrons. 

The measurements are compared to simulations carried out with AliRoot, the ALICE 
event simulation and analysis framework. The non-linearity of the electron drift 
was calculated with GARFIELD and included in the AliRoot code. Also a simpli- 
fied picture of the secondary processes following the transition radiation absorption 
was added to AliRoot. The charge sharing between adjacent pads (pad response 
function) was calculated using an exact method (weighting field formalism). The 
performance of the detector is well understood and the position and angular reso- 
lution are within the requirements for the ALICE TRD. Our results - in particular 
the investigated systematic effects, the corrections applied, and the influence of 
the transition radiation - are of general interest also for other TRD's and/or other 
drift chambers with similar geometry, where a drift region is added to a multiwire 
proportional chamber, with the electron drift perpendicular to the wire planes. 
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